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REVIEW OF DATA ANALYSIS PROCEDURES 
FOR THE ATS-6 MILLIMETER WAVE EXPERIMENT 

Robert Meneghini 

ABSTRACT 

Predictions of satellite downlink attenuation through the use of ground-based measurements 
form a substantial part of the ATS-6 Millimeter Wave Experiment (MWE). At the downlink 
frequencies (20 and 30 GHz), the major causes of attenuation are the density and the size 
distribution of rain drops along the propagation path. Ground station data, which include 
radar and rain gauge records, measure quantities related to the meteorological parameters of 
interest and thereby provide a prediction of downlink attenuation with which the measured 
attenuation can be compared. 

The calibration and data analysis procedures used in the MWE are reviewed with the object 
of improving the accuracy of such ground-based predictions. 

It is hoped that some of the suggested changes can be incorporated into the present data 
analysis procedures. More elaborate data handling and calibration techniques will require 
further evaluation and interim testing to establish their applicability to the Communications 
Technology Satellite and planned reuse of the ATS-6 MWE. 
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REVIEW OF DATA ANALYSIS PROCEDURES 
FOR THE 

ATS-6 MILLIMETER WAVE EXPERIMENT 
CALIBRATION TECHNIQUES 
Introduction 

For the ATS-6 MMW experiment, predictions of the 20 and 30 GHz downlink attenuations 
are found by measurements of rain rate.* The rain rate, in turn, is related to attenuation 
by means of scattering theory as applied to independent spherical scatterers along with ex- 
perimental measurements of rain rate vs drop size distribution. 1,2 ' 

At Rosman, radars at 3.0 and 8.75 GHz, situated near the receiver, are pointed along the 
downlink path (42° with respect to the horizontal). For each radar, an integrated series 
of return pulses (1800) provide an estimate of the return power, and thus reflectivity, at 
locations spaced 100 m apart beginning at a distance 300 m out from the radar. (Note that 
in this report a specific 100 m cell is sometimes designated a range bin or db Z bin). From 
the measured reflectivity (Z), the rain rate (R) can be estimated by a number of experimentally 
or theoretically derived formulas. Details of the radar system can be found in the Westing- 
house reports. 3 

Through the use of ten rain gauge buckets situated under the downlink path, the ground rain 
rate is directly recorded. The data set or its time’ shifted extrapolated version (if the attempt 
is made to account for time and spatial differences between height and ground measurements) 
is then used in the same predictive capacity as that of the radars. 


* An exception is the measurement of sky temperature which, with an appropriate atmospheric model, provides a 
direct prediction of expected attenuation. Radiometric measurements, however, axe not discussed in this report. 



Z-R Relations 


A useful measurement for the analysis of radar data is an absolute calibration of the radar. 
Since radar reflectivity, 2, is used primarily for the determination of rain rate, fixing an • 
empirical relation between these two quantities is essential. The standard relation Z = 200 R 1 - 6 
has certain disadvantages; 


© The variance between the predicted R and the measured R is large, a result due 
not only to averages taken over different systems and different locations but also 
to the inherent variability in the quantity Z being measured. 

© The reflectivity factor Z is a function of the back-scattering cross-section, which 
in turn is a. function of frequency— a dependence which becomes more pronounced 
as the drop sizes increase (i.e. as the deviation from the Rayleigh approximation 
increases). The Z-R relations change with the radar frequency used. In fact using 
the data' of Medhurst 1 , Setzer 2 and Stephen 4 , the-following expressions (derived 
in Appendix A) result 


Z 3 0 =315.5 R 1 * 4 * 7 
Z 87S =307.1 R 1 - 546 


( 1 ) 


• The actual values of the reflectivity factor of 8.75 GHz are obscured by the effects 
of signal attenuation. Again using the Medhurst and Setzer data, the difference 
in. attenuation between 8.75 and 3.0 GHz radar (derived in Appendix B) is 


Ak (db/km) « 0.006569 R 1 - 29 ~ k {ot 

Ak (db/km) « 3.99 X 10* s Z° ®f « k . 
Ak (db/km) * 5 . 544 X 1 O' 5 f ~ k 

0-75. tot 
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The reflectivity factor from the 8.75 GHz radar should be greater than that of the 3.0 GHz 
radar for rain rates greater than about 2 inm/hr. For example, for an R = 50 mm/hr, 
dbZ g75 - dbZ 3i) » 1.4. This statement, however, holds only if the selective signal attenua- 
tion at 8.75 is first taken into account. For example, writing out the radar equations for 
the 8.75 and 3.0 GHz radars: 


8 r 2 




C 3 Z 3 


r 2 


Rearranging 


Z, =r 2 P,/C 3 
and 

Z ' 8 =r2P 8/ C 8 =Z 8 e ' 2/ ° kdr 


(3) 

(4) 


Z 3 and Zg would be the reflectivity factors computed from the radar return power if 

all attenuation effects were ignored. Since the attenuation effect is cumulative it becomes 

more pronounced with an increase in range. At the near range cells (small r) Zg > Z 3 as 

it should be for agreement with rain rate predictions of Equation 1. As range increases, 

-2 f^kdr , 

however, e ° decreases with a corresponding decrease in Z g . For large enough. k, 

Zg < Z 3 beginning at some range r D and continuing thereafter. As a result, Z 3 and Z' will 
lead to differing predictions of rain rate. 

This type of behavior is shown in Figure 1, which plots dbZ vs range. 4 The dbZ' g 
values dominate up to r ~ 1.5 km, after which dbZ' g drops and remains less than dbZ 3 . 


3 




Figure 1 . Storm Cell Structure 

If attenuation effects of the 8.75 GHz signal are taken into account and different Z-R relations 
are used (depending on the radar frequency), the rain rates predicted by the two radars 
exhibit a better agreement than if attenuation effects axe neglected and the same Z-R relation 
is used throughout. 

Figure 2 gives, the rain rates predicted from four-second averages of the radar data. The plot 
marked “Standard Prediction” plots Z = 200 R 1 6 for both radars without attenuation effects. 
The plot marked “Altered Prediction” uses Equations (1) and (2): 

The equations given by Ippolito 6 relate rain rate to downlink attenuation by 

N 

“ 2 o( db ) = £ 0.0687 RR/ 1 (0.1) (5) 

• i=l 
N 

a 30 (db) = 0.1649 RR 1 - 035 (0.1) 

i=l 

N = number of range bins 
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Figure 2. Rain Rate vs Range (Continued) 
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Figure 2. Rain Rates vs Range (Continued) 
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Figure 2. Rain Rates vs Range (Continued) 
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Figure 2. Rain Rates vs Range (Continued) 
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Figure 2. Rain Rates vs Range (Continued) 
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Figure 2, Rain Rates vs Range (Continued) 




Table 1 

Values of Alpha for 0 = 0.887 


GMT 

Attenuation 

Altered 

Prediction 

Standard 

Prediction 

3.0 GHz 
Radar 

8 75 GHz 
Radar , 

ABS. 

Difference 

3.0 GHz 
Radar 

8.75 GHz 
Radar 

ABS. 

Difference 

2323:20 

o(20) 

3 607 

3.525 

082 

3.529 

4.063 

.988 


«(30) 

7.256 

7 082 

.175 

7.113 

8.101 

.534 

2324:00 

o(20) 

3.56 

3 54 

022 

3.56 

3.96 

.40 


a(30) 

7.19 

7.15 

.047 

7.19 

1.95 

.75 

2325.20 

o(20) 

3 287 

3.829 

.542 

3.218 

4.073 

.856 


o(30) 

6.51 

7.48 

.968 

6.338 

7 927 

1.539 

2326-20 

«(20) 

3.178 

3 653 

.475 

3.079 

3.738 

659 


o(30) 

6 229 

7 074 

845 

6 058 

7.225 

1.167 


The Figure 2 plot labeled Altered Prediction (with first 300 m correction) takes into account 
the attenuation introduced by the rain rate in the first 300 meters along the path. The rain 
recorded by the rain gauge buckets has been used to predict the added attenuation along the 
path by 


k(db) = 0.00657 R 1 ' 29 (0.3) 


( 6 ). 


derived in Appendix B. 

For rain rates greater than 1 0 nrm/hr in the first 300 m of path length, the improvement is 
probably significant enough to justify the increase in the necessary computational time. The 
fact that the radars can not measure the first 300 m out and the errors introduced into the pre- 
diction of the downlink attenuation are discussed in the last section of this report. 

Although the above procedures seem to produce some improvement, there are several 
reasons for attempting to calibrate the radars using measured rain rate: 













































<» to confirm or correct the above relations of Z 3 - R, Z g - R 

• to find how the Z-R relations vary with rain type 

• to compare the measured rain rate with the radar predicted rain rate and adjust 
the parameters accordingly. 


Calibration Procedure 


If both radars are pointed in a nearly horizontal direction, x, over the rain gauge buckets 
then the radar equations become: 

p 3 -CaV*’ 




C 8 z 8 C 8 Z g 

8 8 . 8 8 2/k'dx- 


1 0 - 0 . 2 jkix = 


e 


where 


k = 4.343 k' 


Taking 10 X log of both sides, solving for dbZ 3 and dbZ 8 and approximating the integral 

x N 


f k dx by kj, then at range r. 

i=l 


db Z 3 (r.) = dbP 3 (r.) - dbC 3 + 2dbr j 

j 

dbZ 8 (i.) = dbP 8 (r i )-dbC 8 +2dbr j +2 k ( r j) 


(7) 


x=i 


If a measurement of the rain rate is made at rby means of rain gauge bucket then 

R = CZ 3 

r =cz r 

or 
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( 8 ) 


db R = dbc + p dbZ 3 
db R = dbc* +p'dbZ g 


Combining (6) and (7) gives: 


'dbc-dbRR(r.)\ 
dbC 3 =dbP 3 (r j ) + 2dbr j + \ J 


dbC 8 = dbP 8 (ij) 


dbc' - dbRR(r,)\ 
+ 2 dbr j+ ^ j 


(9) 


+ 2 


£ k v 


i=i j 


where k(r.) = k(100) ; + k(200) + + k(j X 100)- ' «. 


i= 1 


and 


k(r) is the one-way attenuation/ 1 00 m at range r. ' 


1 l V 


Equation (9) indicates that once a Z 3 - R relation is assumed, the radar constant C 3 can be 

•f 

found. Alternatively, if C, and p are assumed, c is computed or if C, and c are assumed, 

J * «x w . H * 3 7 

p is computed. In general, for the computation of C g in (9), some estimate of attenuation 
must be used, such as (2); although for small' range, values; and rain rates less than 25 mm/hr, 
the last term on the right side o£(7)*ismegligible and Cg is amenable- to direct computation, 

vfi : ; ! ii’ra 4 i A A V 1 ' 

FUNCTIONAL RELATION BETWEEN SURFACE AND ELEVATED RAIN RATES 


Experimental Procedure 

It was noted previously that rain gauge measurements give a direct record of rain rate on the 
ground and not along the path of interest. Conventional linear or exponential extrapolation 
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from ground upward is usually inadequate because of the temporal and spatial variations of 
■rain. Applicable experimental and theoretical methods are developed herein in order that 
more accurate reconstructions of rain rate at elevated heights can be made. 


The measurement procedure, which is a slight variation of the calibration technique-des- 
cribed previously, employs two radars. One looks along the horizontal over the rain gauge 
buckets' while the other looks along various elevation angles. 

The radar equations can be written as: 


C.Z 


P,(r)=- 


1 1 e -2/k^dr 


•c 1 z 1 


jQ-to/kjdr 




C 2 Z 2 


e-2/k 2 dr 


C 2 Z 2 

-22 io-°- 2 JV r 

r 2 


in which 


k » 4.843 k' 


Assuming that at least one radar is non-attenuating, say radar 2, then: 

W r ) = C 2 Z 2 /r 2 

Next let the non-attenuating radar be directed along an -elevated angle 8 and the attenuating 
radar directed along the horizontal (using the ndtation 'at Figure 3) then: 

r x 

P. (x) = e' 2 Jo k ! dx 




C,Z,( x ,h) 
(x 2 + h 2 ) 
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Writing: 

Z 1 (x,o)=*a 1 R b (x,o) 

Z 2 (x r h) = a 2 R b (x4i) 

where , a 2 can be found either by the calibration technique or some suitable empirical 
formula (Appendix A or B). 


Let 


R' 


'(xh)J =R b (x,o)f(h)J 


where f(h) is the function relating the rain rate at ground (x,o) to the rain rate along the beam 
at the point (x,h). In general, f(h) will be a function of x and h. 



Figure 3. Radar Direction Notation 


Combining the above relations substituting in the equations for return power and taking 
the ratio, then: 


f(h) 




C ! P 2 


C P 


2- X 

— (1 +h 2 /x^)e~ 2 ^o k L dx 
a 2 


An important case occurs when two non-attenuating radars are used in which case k L ~> o 
and: 
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, m 


I 


C P 
^I r 2 

C P 


(H-h 2 /x 2 ) 


Only when the wavelengths of the two radars are identical can a x be set equal to a 2 . This case 
will also, occur if the same radar is -steered from ,the horizontal-to an ; elevation angle d in a time 
which is less than a significant change in the rain rate. Then: 


f(h) 


C P 
12 

C P'- 
J ^2 r l 


C 1 P 2 

2 i„2\ - ^ $ec 20 


(1 +h 2 /x 2 ) = 


C -p 


where 

(1 + h 2 /x 2 ) = sec 2 0 

and 6 is as shown in Figure 3. 

If the above conditions are not met, i.e. if the same A is not used or one is at an attenuating 
A, then a method is needed to evaluate the attenuation factor e" 2 -^ kdx . Possibilities include 

® Use the theoretical expression for K, either (2) or the one given in Appendix B. 

• Use the calibration procedure after first having assumed a Z g 75 - R relation. 

® Consider: 



C 1 a 1 R b (x,o) 


- - e - 2 j> d * 


( 10 ) 
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Taking the derivative with- respect to'range: 


C 1 a 1 R b (x,o) 

dP,/dx= 

1 ^ 


-2 f X kdx 
! J o 


X 


2k L (x) + 

x 


t 

R b (x,o) 


dR b (x,o) 

dx 


(ID 


Dividing (1 1) by (10) and using 


dPj/dx dlnP I 

Pj dx 


and 

then: 


so: 


and 


k(o) = 0 


d InPj 
dx 



-2k(x)+ — lnR b (x,o) 
dx 




C P a 

(1 + h 2 /x 2 ) e- 2k « x o 

C 2 Pia 2 


Thus, f(h) depends on the ground rain rate (through the factor of k(x». This rain rate, 
R(x‘o), can be found by rain gauge buckets positioned at x and x + Ax. 
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These equations follow a procedure similar to that given by Goldhirsh and Katz 10,11,12 . The 
accuracy of such techniques, however, at the frequencies 3.0 and 8.75 GHz do not seem to 
warrant the labor involved. The first two procedures, therefore, are preferable. 

Theoretical Procedure 

A recent personal communication 7 concerning a presentation of the 1975 URSI meeting 8 , 
describes a computational technique for relating ground to height rain rate: 

Let the rain rate be recorded on the ground at times tj , t 2 , . . . t k . The problem is to find the 
rain rate at height h at some previous time t . By means of the Laws and Parsons drop size 
distribution, the percentage of the total water volume contributed within a certain diameter 
increment is known. The number density for a given ground rain rate is given by 

r 

N D =Rp/1.885X 1G 4 v D? 


where: 


Np = density of drops centered about a particular drop size D { . 
R = rain rate 


v = velocity of drops of size D.. 

p = percent of total water volume consisting of drops centered about a 
particular drop size D.. 


Consider the ground rain rate recorded at time t.. Then, at the height h and at the previous 
time of interest, t o , the velocity of those drops centered about a diameter D ; (assuming 
terminal velocity of drops from y = h to y = 0) is given by 

v (D{) = h/(t. - t Q ) 
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If v (Dj) is a known, single-valued function of D. (Medhurst 1 ) then the inverse exists: 

= f(Dj) or D. - f * 1 (*»-) 

Since the rain rate at time t. has been recorded, the density of these drops of diameter 

can be found. Assuming that no change in the drop size structure occured during the drop’s 

downward course, then this is the same density that existed at height h at the previous time, t . 

Similarly for the ground rain rate at another time t , 


v(Dp = h/(t.-t o ) 

and 

^ =f * 1 (^j) 

From the ground rain rate at time t. this densityis found of these drops centered about Dj 
and hence the density at y = h, t = t Q . " 

The entire drop size distribution can be reconstructed at height h and time t G . Assuming that 
the- drops'have reached terminal velocity' at height h then the rain rate can be computed 

•N‘ 

R(y=h,t=t o )= ~ DfN D -(p i ,h)P t (D i ,h). 

i“l 

where N D (D.,h) is the density of drops centered about diameter D. located at the height h 
and i> t (D p h) is the terminal velocity of drops D. at height h. 
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Several criticisms can bemade of this technique. 'It is not.evident that.the terminal velocities 
of drops can be used for drops falling from height h to ground. Furthermore, the expressions 
for terminal velocity assume the absence of updraughts- an assumption which is usually not 
the case especially in severe storms. For cross winds, it is necessary to compensate for the 
horizontal displacement of drops as they fall. 

•Despite these problems, the.technique seems to ;be more rigorous and potentially more 
accurate than the usual time shift and linear extrapolation .used -to -relate ground to elevated 
rain rates. Correctional procedures, such as directional anemometers placed near the rain 
gauge buckets, along with the use of radar as described in “Calibration Techniques” herein 
may -be possible. 

Some preliminary results are given in Figures 4 and 5. The reconstructed rain rates (solid 
unmarked lines) have been superimposed upon Westinghouse curves of rain rate vs, time as 
given by the ground rain buckets (solid curves with crosses) and the appropriate radar bin 
predictions 4 (dotted curves), i.e. that radar bin at a height h directly over the rain gauge 
bucket. The radar data.given in the Figure 4 was taken at a height of 884 m. It should be 
noted that the reconstructed rain rate differs markedly from a time-shifted version of the 
original ground rain rate curve. Figure 5, plotted at 330 m, shows less-pronounced difference; 
a reflection of the fact that total drop size structure for smaller heights approaches the struc- 
ture at ground level. Since the data were taken during a fairly intense storm some measurement 
of updraught and cross wind velocities would be needed for a more accurate reconstruction. 

One other possible procedure for improving accuracy is by the use of the terminal velocities 
of rain drops aloft. 9 It is expected that the accuracy of this reconstruction technique will 
improve as the amount of reliable information input to the model increases. 

TWO-FREQUENCY RADAR TECHNIQUES 

Goldhirsh-Katz Method 

A series of articles* 10 ’ 11 ’ 12 * by Goldhirsh and Katz outline how various storm and radar 
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Figure 4. Rain Rate vs Time at 884 Meters 
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Figure 5; Rain Rate vs Time at 330 Meters 
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parameters can be derived through the use of a multifrequency radar system. In a two fre- 
quency system consisting of an attenuating and non-attenuating radar, the average attenuation 
within a certain range interval is found to be essentially a ratio of radar return powers 
measured at the end points of this interval. 

The equations that follow are the same as those of Goldhirsh and Katz except for elimination 
of the assumption that the effective reflectivity factors for the two radars are identical. The 
only condition necessary is that: 


b=b' 

where 

Zna"*** 

Z A =a'R b ' 


where NA denotes non attenuating 
A denotes attenuating 


Writing the radar return powers of the two radars at the ranges r and r + s 



(12) 

P A (r) = CiZaW 10'»< k4r 

(13) 

r 2 


C 2 ^NA^ r + S ) 

W r + S > = 

(r + s) 2 

(14) 
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( 15 ) 


? 1 ^ A ^ ®) rr + s 

P A (r + s)= — — 10'°- 2 io kdr 

(r + s ) 2 

Dividing ,(12) by (13) taking 10 X log of the resultant and using -the relations Z NA = a 2 R b 
and Z A = R b 


then: 



P maCO 



f c z\ M r 

10 log 


— 1.0 

logj 

— ) flog H- +0.2r kdr 

\cj AV ' J 


(16) 


Proceeding similarly with (14) and (15) 


10 log 


?NA< r + S > 

P A (r+‘s) 




(17) 


Subtracting (16) from (17) and approximating 
average, then: 



"* r 


fey ks where the bar denotes an 


1 

k = ~ 

2s 


** bI W r + ?)* #?na( 1 > 1 4 bP A (r) r dbP A (r f. s) 

- 


(18) 


where s is expressed in km and k is the average attenuation in the interval s in dfe/km. 


The advantages of equation (18) are the independence of k" on the radar calibration constant 
and the absolute range, (i.e. only the increment length s enters in). The results are due to. 
the fact that k is sensitive only to retwn power changes in the Interval (r, r+ s). One other 
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advantage of the method is that- if the exponents in the Z-R relation are taken to be equal, 

(b = b f ), then K is also independent of the coefficients \ and a 2 in the Z-R relations. 

The drawback of this and similar techniques are the errors arising from the statistical nature 
of the target. The fluctuations are caused by changes in the relative positions of the particles 
(i.e. positional changes of a significant part of the radar wavelength) and result in variations 
of the incoherent addition of backscattered powers from the individual scatterers. Such 
changes occur in times on the order of milliseconds. 

Among the parameters necessary for the determination of the error variance of Z are the sig- 
nal to noise ratio and the number of dependent samples used to estimate Z 14 . It should be 
noticed that for a given Z value the errors will depend upon range and attenuation as well as 
the magnitude of Z as these quantities determine the backscattered signal level. 

For an estimate of the number of independent samples taken by the radars, an equation from 
Battan is used 

r 01 =1.71 X lO' 3 sec 

where t 01 is the time required for the autocorrelation of the signal to fall to 01. 

For the Rosman radars, the averaging is done over 1.8 sec at a PRF = 1 KHz. Therefore, for 

#• 

the 3 & 8.75 GHz radars 105 and 307 independent samples, respectively, are obtained per 
range bin out of a total of 1800 return pulses. Noise level inputs to the preamp are -102 dbm 
and -99.5 dbm for the 3 and 8.75 GHz radars respectively, while the signal levels are found 
from the standard radar equations. 

Statistical analyses required to complete error variances of reflectivity and rain rate or to 
assess the viability of multi-frequency radar techniques could not be included in this report. 
Nevertheless, it is suggested that statistical analyses be taken as an essential part of future 
radar data processing. 
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The method used to analyze the ATS-6 data using the Golhirsh-Katz method is given below. 


For a given increment length S, equation (1 8) is used to compute attenuation, using an 
average of three data points as estimates of return power at the ends of the interval. Next, 
the average dbZ 3 , value within the increment S is computed. As a result, pairs of values (k, 
dbZ) are found throughout the entire data set of radar return power vs range. 

The final step is finding the average and variance of those values of k which correspond to 
dbZ values falling in some AdbZ increment centered about a particular dbZ. The entire process 
is repeated for various incremental lengths S ranging from 0.5 to 1.3km. Some of the results 
are contained in Table 2 below in tabular form where the values of k as computed by the 
Goldhirsh-Katz method are compared with values of k (in db/km) computed by the method 
of Appendix B. 


Table 2 

Comparative Attenuation Results 


DAY: 270 GMT. 2325:20 S=.9km 


dbz av 

k 

(Goldhirsh-Katz) 

k 

= 0 00029 Z- 72 

k 

= 0 00004 Z* 51 

39.9 

0 156 

0.22 

0141 

40 4 

0 177 

0.235 

0.153 

40.7 

0 138 

0 247 

0.163 

41.3 

0.433 

0.273 

0.184 

41.5 

0.242 

0.283 

0.193 

42.6 

0997 

0.339 

0.24 

42.8 

0.285 

0 349 

0 249 

43.25 

1.12 

0.377 

0 247 

43.7 

0.339 

0.406 

0 30 

44.1 

0 368 

0,432 

0 325 


General characteristics of the results are: 


• As S increases, for a given number of points averaged, the variance of k decreases. 


® Although the results show that an increase in dbZ is usually accompanied by an increase 
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in k, for a significant portion of the data k is not a monotonicaUy increasing function 
of dbZ. 

© For the data analyzed here the Goldhirsh-Katz method usually predicts a larger 
value of k than the theoretical formulas. 

In the absence of an adequate error analysis, the accuracy of the Goldhirsh-Katz technique is 
difficult to assess. In general, the predictions are larger than those given by the theoretical 
formulas. Figure 1 demonstrates that if the 8.75 GHz calibration constant is decreased, the 
8.75 GHz signal attenuation must increase for proper agreement of predicted rain rates. That 
the 8.75 GHz calibration constant may be too large can be seen from the curves at Figure 6 
which compare ground rain rate with the corresponding 8.75 GHz predictions. The height 
difference between rain bucket and range bin is fairly small, i.e. 102 m so a correlation co- 
efficient near 1 is expected. It should be kept in mind, however, that for the computation 
of the rain rate, the same Z-R relation (Z = 200 R 1 - 6 ) was used throughout. 

In addition to experimental measurements of attenuation, Goldhirsh and Katz demonstrate 
that, in principle, an estimate of drop size distribution can be found by employing an attenu- 
ating and non-attenuating radar. Once this distribution is known over the path of the down- 
link signal, the expected attenuation of the 20 and 30 GHz downlinks can be computed with- 
out having to deal with the intermediate quantity of rain rate. 

Assuming the Marshall-Palmer form for the drop size distribution, 

n(D) dD = n o e' AD dD 

and 

Zna-JT 5 NA 0»n(D)D 6 dD 

and 

k=/:C ext (D)n(D)dD 


(19) 

( 20 ) 
( 21 ) 
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TIME IN MINUTES 


START TIME STOP TIME 


YEAR = 74 YEAR = 74 

DAY = 270 DAY = 270 

GMT = 1910 GMT = 1940 

NO SECS AVGD = 20 


RADAR BIN - 1 
RAIN BUCKET = 2 


Figure 6a. 3.0 GHz Radar vs Rainbucket Rain Rate 
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RAIN RATE (mm/hr) 



TIME IN MINUTES 


START TIME STOP TIME 


YEAR = 74 YEAR 

DAY = 270 DAY 

GMT = 1910 GMT 

NO SECS AVGD = 20 

RADAR BIN = 1 

RAIN BUCKET = 2 


74 

270 

1940 


Figure 6b. 8.75 GHz Radar vs Rainbucket Rain Rate 
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•and taking the ratio k/Z NA and using equations (19), (20), and (21) then 


k f: C.„,(DK AD dD 
Z NA jf*o NA (D) e' AD D‘dD 


where Z NA is given by equation (12), k by equation (1 8) and -values of C ext (D), o NA (D) 
obtained from numerical computations using Mie scatter theory. 1 ’ 2 

Equation (22) can be solved numerically for A which, in:tum,is-used to, compute- n o 
through equation (20) 


a. 


= Zna//o* e' AD D 6 dD 


(23) 


In a.recent report on error analysis- of- this technique 13 , an,analytic,.expression ismsed to 
approximate (22) and- (23)- for two radars at 3 and flO'GHz. Itis;found that: 


and 


A(cm* f) = exp 1 33 .094 - 54.267 x 7.4182 x 2 - 0.-3322 x 3 | 

l i 


(24) 


n o =7 J na/ cx 1 ) ! *2.95 -6.16 In A- - 0:lil'32(ln A) 2 j 


*(25) 


where 



From Battan 5 , pg 44, Z 0 , the effective reflectivity factor, is related to - tj,! the. reflectivity, 
through the equation: 

Z e =X 4 ij/ir s :KP 
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where 


m 2 - 1 


«« 0.93 


lkP = 

m 2 +2 

m = complex index of refraction for water 
X = radar wavelength 


for 

X 3 = 10 cm 

t? na (km' 1 ) - Z 3 X (2.846 X 10’ 9 ) 


where Z 3 is in units of mm 6 /m 3 


x 



If a Z-R relation is assumed, such as the one given in Appendix A, 


R = (Z 3 /261) 2 ^ 


then N d vs D can be plotted for various values of rain rate and attenuation, k. Since the 
Goldhirsh-Katz method provides neither a convenient analytic representation of k nor re- 
liable values of it for radar frequencies of interest unless a considerable amount of averaging 
is used, an alternative approach to finding attenuation is considered in the following section. 

Alternative Technique 


For the experimental determination of radar attentuation, a procedure somewhat different 
from the Goldhirsh-Katz technique is considered. The motivation is that the predicted rain 
rate from each of the two radars should be the same and that the discrepancy between the 
predicted values should be assigned to the selective attenuation of the higher frequency radar. 
Equating rain rates, or equivalently, 10 log R, 

10 log R 3 ' = 10 log R 8 

10 log(aZ 3 ) = 10 log (a'Zjj*) 
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or 


10 log 



■H0*b*logZ 3 = 10 • b' • logZ 8 


If b = b', then, at a single range cell 
where 


K + dbZ,.=dbZ n 
3j j 

K = 101og^/b 


( 26 ) 


The actual value of dbZ g is the sum of the ‘measured’ dbZ g plus the amount of attenuation 
suffered through the two-way path, where the ‘measured’ and actual values of dbZ g are 
defined by 


CZ 8 -r , ^8m 

p= — i — io' 0-2 4 k ^ = — — 
r 2 r 2 


in which the subscript ‘m’ denotes measured. Taking 10 X log of each of these equations 
and Subtracting one from the other 


d « 5 .=i!=Z s „. + 2 £ k i 

i=l 


(27) 


where 



Relating k. to Z 3 by the equation 

!c =«Z» 

equation'(27) becomes: 

i 

dbz s j =dbz 8 »i +2 Y 

i=i 
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Equation (26) can be rewritten as: 


j 

K+dbZ 3j =dbZ ln , J+ 2 J2 12* 

i=l 

This equation characterizes one range bin. For N range bins, N equations result which, when 
added, give: 

N 

‘ X dM »mJ +2 « 

j=*'l J“.l i=l 

* ‘ a 1 

The easiest procedure is to fix 0 and solve for a from which: 

N 

NK+ X ( db V dbZ W 

«- FT 7 (28) 

X X z ?, 

j=l i=l 

Note that for computational ease the denominator can be written 

N j 

X/ Z f =NZ 3 1 + ( N - 1 ) Z 3 2 + f N - 2 ) Z f 3 + --- +Z 3 n 

j=l i=l 

N 

- X (N - bti ) z ? k 

k=l 

The attenuation becomes 

k = aZ <5 (29) 


X 
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where a is given by equation (28). 


Equation (29) can be interpreted as the attenuation of the 8.75 GHz signal such that the 
average absolute difference in the predicted rain rates is minimized. Unlike the Goldhirsh-Katz 
equation, which computes attenuation in an interval, equations (28) and (29) characterize 
the set of data points as a whole. Another difference between the two methods is that a, 
and hence k, depend upon the Z-R relation; a relation which depends upon the radar calibra- 
tion constants. Nevertheless, this technique is an easy means of comparing the theoretical 
expressions for attenuation with an expression found from the data. If such a procedure 
is carried out over a sufficient amount of radar data, characteristic values of a are found 
for a number of different classes of rain. 


Computed values of a are tabulated in Table 3, in which 

k - ctZ 0 " 88 

These values of oc should be compared with the theoretical expressions 

ic = 0.00004 z 0 - 88 


and 


k = 0.00029 Z 0 ' 72 

The sharp decrease in alpha between times 2324: 10 and 2324:30 corresponds to a fairly 
rapid movement of the storm towards the receiver, a general increase in rain intensity in 
the near range cells, and a decrease in the spatial extent of the storm. The increase 
around GMT = 2328 corresponds to a decrease in rain intensity with the spatial extent of 
the storm basically unchanged. From time 2324:30 to 2327:50, the computed values ■ 
are in fair agreement with the theoretical value of 4 X 10 s . For the remaining intervals 
2323:10 to 2324: 10 and 2328:10 to 2329:50, the value are greater than two and in 
some cases as large as three times that given by theory. 
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Table 3 

Phase Changes in Different Rain Rates 


GMT 

ALPHA 

GMT 

ALPHA 

2323:10 

.000129 

2326:10- 

.0000634 

2323:30 

.0001 13 

2326:30 

.0000408 

2323:50 

.000123 

2326:50 

.0000457 

2324:10 

.0000989 

2327:10 

.0000336 

2324:30 

.0000523 

2327:30 

.0000535 



2327:50 

.0000635 

2324:50 

.0000431 

2328:10 

.000107 

2325:10 

.0000498 

2328:30 

.0000824 

2325:30 

.0000355 

2328:50 

.0000977 



2328:50 

.0000977 

2325:50 

.0000486 

2329:30 

.0000939 



2329:50 

.000129 


One should not conclude, however, that the theoretical expressions for attenuation are 
necessarily in error by taking as evidence the results of both the Goldhirsh-Katz and the present 
method. In fact, the Goldhirsh-Katz technique often gives highly variable and, therefore, sus- 
pect results while the present method might he using incorrect values for the radar calibra- 
tion constants and therefore also be in error. Even if it is clear that the theoretical formulas 
underpredict in certain cases, the problem of explaining the mechanism behind such changes in alpha 
remains unsolved. One possible cause for the variation in alpha may be a shift in the drop size 
distribution toward larger drop sizes thereby accounting for an increase in attenuation and 
therefore in alpha. Disdrometer data, along with radar calibration curves should be helpful. 

It should be noted that any anamoly in the attenuation data for the 8.75 GHz radar most likely 
has an analogue in the measured attenuation of the 20 and 30 GHz downlinks. Therefore a pro- 
per explanation of the radar attenuation should also help to resolve the discrepancies between 
the theoretical and measured- downlink attenuation. 
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Having computed several values of a, theoretical drop size distributions can be plotted based 
on equations (24) and (25). and the Z-R relation Z 3 = 261 RR 1S . Plots of N D vs D for 
various values of rain rate and attenuation are presented in Figure 7. Also plotted on each 
graph is the Marshall-Palmer distribution (solid lines) given by: 

Njj(m' 3 mm' 1 ) = N o e" AD 

where 

A (cm' 1 ) = 41 RR' 0 21 

N 0 (cm' 4 ) = .08 
D - drop diameter 

The rain rates that correspond to the solid lines are 1, 5, 25 and 50 reading up in a counter- 
clockwise direction. Certain characteristics of the curves seem to be qualitatively correct. 

For example: (1) increasing attenuation is accompanied by a greater percentage of larger 
drops/ This accounts for attenuation increase because for a given volume of water, the volume 
composed of larger drops results in more attenuation; (2) N D decreases exponentially as D 
increases, and' (3) for small D, N D increases for larger fain rates, which is in agreement with 
the modified-Marshall-Palmer distribution. 5 

A disturbing feature of the curves is that for a > 0.000042, the equations predict a greater 
percentage of larger drops for rain of 25 mm/hr than at 50 mm/hr. One possible reason is 
that «> 4 X 10' 5 is outside the region of validity for the expressions used. For a < 4 X 1 0“ s 
the curves appear qualitatively plausible and it might be of interest to compare height corrected 
disdrometer data with curves such as those given here. 

Due to calibration errors and the statistical fluctuations of scattering from hydrometers, 
accurate drop size distributions from radar measurements probably cannot be obtained 
without going to a higher .frequency radar. Nevertheless, experiments or error analyses would 
be of interest not only because of the theoretical importance of such techniques but also for 
assessing the feasibility of a higher frequency radar for future systems. , 
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Figure 7. Rain Rate vs Drop Size (Continued) 
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Figure 7. Rain Rate vs Drop Size' (Continued) 
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PHASE DIFFERENCE INFORMATION RELATED TO METEOROLOGICAL PARAMETERS 


At Rosman, differential phase measurements are made between the carrier and various side 
tones: The question of interest is how this quantity varies with the meteorological conditions 
along the downlink path. The procedure, following Van de Hulst 15 , uses an effective refraction 
index to characterize forward scattering through a multi-particle region. 


Van de Hulst considers a slab of length £ filled with N particles/vol, identical and identically 
oriented. The distance Z from the Observation point to the slab is assumed to be much larger 
than a wavelength, X. Further, if the effective refractive index is such that m « 1 , then the 
ratio of total to incident disturbance is: 


_ _ e -ik£(m-i) _ 


2 7T 

1 NVS(O) 

k 2 


let 


m = n - in' 


then 


- 1) = 1 - ilc2( n - in' - i) 


Therefore 


n = 1 -r 2?rNk" 3 Im 


S(0) 


(30) 


where S(o) denotes the amplitude scattering factor in the forward direction. 


Equation (30) determines the phase lag of the wave. The effect of such a phase lag, as 
Van de Hulst remarks, is a dispersive medium characterized by a wave velocity v, where 
v - c/n. 


If the particles have a distribution of radii given by N(a),- then: 

~ > i2ir 

m = n - in =1 /S(o, a) N (a) da 

k 3 
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and 


n = 1 + — /“ Im S(o,a) | N(z,a) da 


The drop size dependence on z is used simply to emphasize the spatial variation of the dis- 
tribution. The quantity S(o,a) can be related to the Mie scatter coefficients by 


so that: 


oo 

S(o,a)= — - £ (2n+l)(a n +b n ) 


n=l + — £ S (2n+l)Im(a n +b n )N(z 3 a)da 



Figure 8. Notation for Path Length 


Referring to Figure 8 and using the relation 


n = cjv = c/df/dt 


t= — £ na)df 
c 


Substituting equation (31) into (32): 


oo 

+ - £ C E (2n+l)Im(a n +b n )N(z 3 a)dad? 


C 
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The phase of the wave can be written: 

k' {:. s; ’ \ f- '■"/ 

\ 

(,s - 2i:(t/T); T = period 

•r-> ' t J«i n -t y !' i- .t’-’i’b 1 

For the ATS- 6 Millimeter Wave experiment, three monochromatic .signals are sentdn phase 

,u •.,)!* "■ tv' '.>.>'*• < 7-« r • i ■ 

from the spacecraft. The quantity measured at the ground station is: 

4 ' i 


» i 


1 . 


Where the subscripts C, U, L denote the carrier and the upper and lower sidebands. 
For clear sky measurements, the 'above quantity is calibrated to zerp: ‘ " 


y’dNR ~ 





( 34 ) 


Using .equation (34), the differential phase measured in rain can b.e written: 



U 



-0 



(35) 


- ,i' M ‘ 1 . '< ^*.4.; } • < (■>#■». 4 

which is just the phase lag measured with respect to clear sky conditions. Using equation 
(32) gives: 




^dR 


2tt 2 

£ s; E 

n = i 

V k O 

tt °u 

U 1 


•Nti. 

2a 2 

AC E 

n-l 

T» ck 2 

L - o L 

An 2 , 

t . *■' f 

op 

J o wo Z ^ 

n=l 

T ck 2 

c °o 


n ( a nL +b nL )Ndadf &6) 

,{ ! * i * *?m t i i’ 4 t *■ >' hi » - '* 
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Using k = w/c and w = 2tt/T and reducing 


^dR 


4 IT 


2 pntl))< Im K" T » +a »t T i- 2a ~ T =) 07 ) 


n=l 


+ 0i„j;+b nL T J L -2b ne Tjj] XN(z,a)dadS; 


The coefficients a n , b n for the frequencies that are of interest here have not been tabulated. 
A good approximation can be made, however, if the a n ’s and b n ’s are expanded in a Taylor 


series about the center frequency f c 

da n 

a n(0-a n (f c )* 1T 

1 d2a n 

(f-f c )+~ 

2 df 2 

f=f U1 

C 

, db n 

b a( f ) = b „( f c) + 

1 d2b n 

2 df 2 




f=f 




f=f 


C C 

Substituting the first two terms of these expansions into equation (37) gives: 


L cnbi) 


X 


n=l 


fT 2 + T 2 1 2T 2 ) (a a + b n ) 


^ da ni db nn 

+ 1 — - + - 
df df 


[(f„-f t )T; + (f t -f c )T 2 L ] 


(38) 


XN(z^)dadf 


where 


a ni = Im a. b . = Im b „ 

m nc 3 m HO 


Rewriting equation (38) by noticing that only the coefficients a n , b n are frequency-dependent 
gives: 
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*-«- Z C Zl (2n+l)(a m +b » )<T « * T ! - 2T c) 

X N(Z,a)da df 

+ jf 4 C 2p»+1Kv + '>J [(C-fX + < f L-WTl] (39) 

X N(z, a), da df 

Some idea of the magnitude of fdr as a function of rain rate can be found from the results in a 
paper by Setzer 2 in which a numerical solution is carried out for expressions of the form: 

c 2 N 

^ — So - C X (2n+l)(a ni + b ni )N(z,a)dad? 

4ir ~ 

n~I 

where T k = period. 

By fixing the rain rate and varying the frequency, Setzer’s data show that the second term of 
equation (39) is much less than the first. The sole difference between equations (39) and 
(40) is in the multiplicative constants T 2 and (T 2 + T 2 - 2T 2 ). Thus given the above 
approximations, Setzer’s data can be used to compute y?dr of equation (39). 

For the 20 GHz carrier with ±720 MHz sidebands, 

(T 2 + T 2 - 2T 2 ) = 19.52 X 10' 24 sec 2 

U 1# c 

Setzer compiles data for 

T 2 ~ 2.92 X l<r 21 sec 2 

From Setzer 

T 2 X = v j ( rr > 
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but 


C r2 +T 2_ 2T 2)X= ¥) (RR) 


so 


( T 2 +T 2 _ 2T 2) 

^ dR (RR) v<RR) 

T c 

«p dR (RR) = 0-0067 <p(RR) 


(41) 


where: 


RR = Rain rate 


X= — W y'(2n+l)(a ni + b ni )N(z,a)dadg 
47T 


i p (RR) = Phase change in deg/km 


tabulated by Setzer as a function of rain rate. 


For the 30 GHz carrier with ±720 MHz sidebands, a similar procedure gives: 


<£ dR (RR) = °.0°35 tfRR) 


(42) 


Table 4 lists the values found from equations (41) and (42) for various rain rates. 
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Table 4 

Rain Rate Values. 


Rain Rate (mm/hr) 

<p dR deg/km 

f c -20GHz - 

f = 30 GHz 

.25 

0 

0 

1.25 

.0147 

•0124 

2.5 

.044 

.025 

5.0 

.074 

.017 

12.5 

.148 

.0996 

25.0 

.253 

.175 

50.0 

.431 

.262 

100.0 

.111 

,45 

. 150.0 

1.1 

,612 


It appears that even for extremely large rain rates over a large distance, (say 4 km), the 
total differential phase shifts predicted are: 

¥> dR (f. - 20 GHz, RR = 150 mm/hr) =» 4.4 ? 
yJ dR (f, = 30 GHz, RR =150 mm/hr) =» 2.5° 

The values in Table 4 should in no sense be presumed accurate. Within the accuracies of the 
model, however, a few conclusions might be drawn: 

® The values of the measured differential phase are expected to be small, 

9 <4^ (f, = 20 GHz) > ^ (f c =30 GHz) for all rain rates of interest. 
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• If the actual differential phase is in general agreement with that predicted then 
almost all relevant information will be lost since the system is not expected to be 
sensitive to changes in phase on the order of a degree or less.' 

Some differential phase data have been obtained at Rosman, N.C. on day 270, 1974 during 
a fairly intense storm. The data have been obtained' from Reference 4 and are. displayed 
in Figure 9. Unfortunately, the data seem inconclusive since, over the relevant portions- of 

i 

time, the phase varies within approximately ±4° of a constant value and presumably within 
the errors of the system. The strong phase fluctuations near 2324 and 2328 should probably 
be discounted as the receiver was just losing and regaining phase lock at these times. 

DISCUSSION 

If a few changes are made in the data analysis procedure, the rain rates predicted from the 
3.0 and 8.75 GHz radars should be in better agreement. As a basis for absolute comparison 
between predicted and measured rain rates, the calibration procedure may be helpful as could 
the measured and theoretical relations of ground-to-elevated rain rates. 

Several other changes in the data analysis might be helpful: 

• Range bins 4, 13, 14, 15 for the 3 GHz radar are in error by about 9.5 db. For these 

points, and points such that the 8.75 GHz data exists and not the 3.0 GHz data, 
the later can be reconstructed by use of the equations in Appendix A. The ( - 
estimate of,dbZ 3 at the m + 1 range bin is: * ‘ 

f 

dbZ 3 (m + 1) = dbZ 8 (m + 1) -dbZ g (m) + dbZ 3 (m) + K X .2 t 

where k X .2 is the attenuation suffered by the 8.75 GHz radar between range bins 
m and m + 1 . Normally this last term is negligible. 
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Figure 9. 720 MHz Sideband Phase vs Time (Continued) 



230 _24p;240 250 260 ,270 ;280 290 .300 '310 320 

TIME IN SECONDS 

;2324 '2325 


330 340 350 360 360 


. START-TIME ; STOPTIIVIE 
f.OyTER INTERVAL 


YEAR 

rDAY 

"GMT 

= <1974 
“ J270 
= 2315 

•YEAR 

fDAY 

GMT 

= "1 974 
.= .271 

j= -20 

-TEST MODE 

,20 GHZ 30 GHZ 
MT MT 

DAY 

INNER INTERVAL 
= 270 DAY 

= ;270 

PAR AMP 
-$/C'ANT» GAIN 

OFF 

HI 

OFF 

HI 

,GMT 

= ,2320 

vGMT 

>= 2349 * 

MULT. CHAIN 

A 

A 


Figure 9. 720 MHz Sideband Phase vs Time (Continued) 
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Figure 9. 720 MHz Sideband Phase vs Time (Continued) 
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Figure 9. 720 MHz Sideband Phase vs Time (Continued) 



The radars cannot measure rain rates of distances less than 300 m from the 
receiver. Nevertheless, an account should be made of the first 300 m either by rain 
bucket readings or by an extrapolation from the first dbZ bin. Some idea of errors 
introduced if these points are neglected are 


G 30 

~ 1.7 db 

R = 30 mm/hr 

a 30 

~ 2.25 db 

R = 40 mm/hr 

a 20 

~ .87 db 

R = 30 mm/hr 

a 20 

~ 1.2 db 

R =40 mm/hr 


For those bins (at a given time) with no radar data, which occur between bins 
with data, lower and upper bounds for rain rates in these intervals can be assigned. 
Evidently, the lower bound is zero while the upper bound is the rain rate which 
corresponds to the minimum detectable reflectivity at that range. Thus, in the 
prediction of 20 and 30 GHz downlink attenuations, upper and lower bound 
predictions would occur. Note that in most instances the predictions will not 
differ by much more than 1 db. 

Due to the rain drop’s deviation from sphericity, there is an added effective 
attenuation of the downlinks not accounted for by the radars. To correct for 
this effect, this attenuation due to depolarization must-be added to the radar 
predicted attenuations. The difference, however, is slight. 

Although drop size distribution measurements do not at this time seem possible 
using the Goldhirsh-Katz or similar techniques, the alternative method Two 
Frequency Radar Technique seems to provide a fairly simple means of determining 
the average attenuation of the higher frequency radar. 



• In the millimeter wave literature, no account has been found which describes the 
effects of scatterers in the vicinity of the receiver. It would be of interest to deter- 
mine experimentally whether a positive correlation exists between high intensity 
rain rates in the near fields of the downlink receiver and errors between predicted 
and measured attenuation values. 

» Error analyses of Z due to target fluctuations could be incorporated into the 

present computer programs. Occasional measurements of noise' and clutter levels 
are necessary inputs to the statistical model. Such procedures would impose limits 
on the maximum - possible accuracy of radar predictions of downlink attenuations. 
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APPENDIX A 
Z-R RELATIONS 


Y] N D. D f 5 3. / N(D) D 6 5 3 (D)dD mm 6 /m 3 
i 

T*! Np Dfg^ lim / N(D) D 6 5 g (D) dD mm 6 /m 3 

where „ ^mie 

5= 

^rayl 

is the ratio of the backscattering cross-sections as predicted by the Mie and Rayleigh 
theories. Note that 5 is a function not only of X and D but also of temperature. The cal- 
culations use a rain temperature of 1 8°C. The values of § 3 , S 8 are taken from Reference 3 

I 1 

of the main report. Velocity vs. diameter and rain rate vs diameter are taken from 
Reference 2. 

The expression for N o is given in Reference 1 of the main report. 



N D (in 3 ) - 


RR X (% tot vol) 
1.885 vD 3 


in which 


RR 

= mm/hr 

V 

= m/sec 

D 

= cm 


Computing Z 3j0 and Z g7S for various rain rates and fitting this data to an aR b form by a 
least squares fit then 


63 



Z 3 .0 4 315.5 R 1 ' 457 
Z 8 7S = 307.1 R- 1 - 546 


If b is fixed to 1 .5 then the best fit becomes: 


Z 30 = 261.16 R 1 * 5 
Z 8 _ 75 = 341.715 R 1 * 5 


These values should be compared to those given in Reference 5 of the main report. 


Marshall - Palmer (0°C) 

\ M s - 275 R 1 - 55 


Z 3.0 = 2!0R“ 


Gunn & East (18°C) 


Z 9.34S - 310R 1S « 


Z,.o - 210 R 1 ' 6 


Mueller & Jones (0°C) 


■“ 9 . 345 ' 


3.0 


890 R 
810 R 
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APPENDIX B 
K - Z. RELATIONS 


Ak (db/km) can be found directly from Medhurst’s data. A mean square fit is then performed: 

Ak (db/km) = .006569 R 129 . 

Ak (db/km) » 3.99 X 10* 5 Z°- 88 
Ak (db/km) = 5.5444 X 1 O' 5 Z® 88 

where 

Ak (db/km) = att (Z g „) db - att (Z 3 0 ) db ~ att (Z g7s ) db 

The plots: 


Ak 

vs R 

Ak 

VS Z 3.0 

7 

0.75 

vs R 

Z 3.0 

vs R 


are given on the following pages. 
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K (dB/km) 







K (dB/km) 



Z (3.0) 








